The hyperfine transition in light muonic atoms of odd Z 
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The hyperfine (hf) transition rates for muonic atoms have been re-measured for select light 
nuclei, using neutron detectors to evaluate the time dependence of muon capture. For ^^F = 
, 5.6 (2) fj,s~^ for the hf transition rate, a value which is considerably more accurate than previous 

• measurements. Results are also reported for Na, Al, P, CI, and K; that result for P is the first 
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positive identification. 
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I. INTRODUCTION 



There is still much uncertainty about the hyperfine transition rate for muonic atoms in the Is ground state. Differ- 
ent techniques often give different rates, and there are many inconsistencies. This topic was boosted into prominence 
' 30 years ago by the Chicago group of Winston, Telegdi and co-workers who studied ^^F and established the 
. major properties of this effect. Even today, ^^F remains the best understood example, because of the convenient time 
C3 ' constant of about 180 ns. 



^ When a /i^ stops in a target, it forms a muonic atom, and quickly cascades down to the Is level on a time scale 

^ [ of 10~^^ s which is effectively instantaneous for any detector. If the nucleus has no spin, there is a single ground 

D I state and the muon awaits its fate via the weak interaction decay (/i^ — > e~ Fe z/^ ) or via nuclear capture. For 

-H ' the light elements this occurs within a few /xs and can be studied with a variety of detectors. If the nucleus has 

^ , a non-zero spin however, the situation is more complicated because there are two hyperfine levels of the Is state, 



separated by an energy varying between a few eV and a keV or so. If the transition between these levels occurred 
via an Ml photon, the rate is too slow to be observed. However, the Chicago group showed that Auger emission can 
^ ' speed up the transition rate a thousand fold and bring it to the time scale of nanoseconds, which is within the range 
k>( [ of standard detectors. The hyperfine transition can be detected in any nucleus by observing the depolarization of 
. the /i~ via the detection of its decay electron. However, the ^~ has a very small residual polarization in the atomic 
state (< 10% for nuclei with spin), and in addition the fi~ can be depolarized by magnetic interactions or by nearby 
radicals. Nevertheless, several experiments have been carried out successfully. 

If the nucleus has an odd Z, the is magnetically coupled to the odd proton. Now the capture probability in the 
//~p system, is 660 for the singlet state, but only 12 s~^ for the triplet state. Thus in one hyperfine state the total 
muon capture probability is approximately proportional to (1 -I- Z/2) and in the other proportional to (Z/2). Thus 
nuclear capture is highly sensitive to the hyperfine state, and the hyperfine transition can be followed by detecting 
neutrons from the capture events. The effect can also be observed via the time-dependence of the decay electron, but 
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in light nuclei the normal decay rate dominates and thus dilutes the signal. Winston |2| as well as Suzuki et al. g] 
were able to observe this effect in "'^^F, but for other nuclei the signal is too small. 

The hyperfine transition can also be observed by detecting specific 7-rays resulting from nuclear capture. Some 
transitions are highly sensitive to the initial hyperfine state and, for certain spin combinations, this is also true for 
even-Z nuclei such as ^^C. Several nuclei have been studied by Gorringe and coworkers In addition Wiaux j|] 

has found a very large sensitivity for the 320 keV 7-ray from muon capture in ^^B. An unexpected observation by 
Gorringe et al., was that the hyperfine rate in metallic sodium was Ah = 15.5 ± 1.1 ^s~^ but when in NaF, 
the sodium hyperfine rate was only Ah = 8.4 ± 1.9 ^s~^ the fluorine hyperfine rate was Ah = 4.9 ± 1.2 /is^^, 
consistent with Winston's measurement of Ah = 6.1 ± 0.7 /is~^ ||], and our own measurement of Ah = 5.6 ± 0.2 
lis~^ in LiF. This raised the possibility that there might be a difference between metals and insulators, or between 
different molecular species. Such a difference could be due to different recombination times, or an unusual sensitivity 
to the binding energy of the valence electrons. However both possibilities seem unlikely. 

The most puzzling result to date has been the observation of a hyperfine transition in ^^N by Ishida et al. , using 
the /i~SR technique to observe the muon depolarization. They detected a relaxation rate of 0.092(33) /is~^ of which 
about 0.016 /is~^ is due to the different loss rates from the hyperfine states. This leaves 0.076(33) /is~^ which may 
be from a hyperfine transition. The problem is that the energy difference between the hyperfine states is only 7.4 eV, 
whereas, in the carbon atom, the least bound electron is bound by 11.3 eV, i.e. this energy is the ionization potential 
(the fi~ is close to the nucleus, so the ^'^N muonic system appears to the electrons to be more like a ^''C nucleus). If 
the pseudo ^''C atom has formed a "C"N bond, the ionization potential is 14.3 eV, which makes the matter worse. 
We note that Wiaux found a slower hyperfine rate in -'^^B in comparision to a /iSR experiment, so we suggest that the 
observation in ^"^N was caused by an additional depolarization mechanism, not the hyperfine transition. 

Because of the various inconsistencies in this field, it was decided to re-investigate several light nuclei, using the 
more reliable technique of neutron detection. A comparison was also desirable between Na and NaH, Al and LiAlIl4, 
and K and KH, to see if any difference could be detected. Unfortunately, no satisfactory result was obtained for KH. 
It should be noted that any muon captured in hydrogen is immediately transferred to another atom because the /ip 
system is neutral. Thus NaH, for example, is effecively pure sodium in an insulating environment. 



II. EXPERIMENTAL TECHNIQUE 



The experiment was carried out in the M9B channel at TRIUMF, which includes a superconducting solenoid. For 
muons of 60 MeV/c the stop rate was ~ 10"* with an electron contamination of 2 % and a pion contamination 
of < 0.2 %. The muons passed through 2 plastic scintillators, and stopped in various targets; a third large plastic 
scintillator was used as a veto to define a muon stopped in the target. A mu-metal shield around the targets reduced 
the ambient magnetic field from 1.5 to 0.1 gauss. Table | lists the properties of the targets. 

The neutron detectors were four cylindrical liquid scintillators, 2 of NE213, 1 of NE224 and 1 of BC501A (equivalent 
to NE213). They were arranged in a symmetrical array at 45°, 135°, 225°, and 315° to the beam, in order to minimize 
muon rotation effects. Plastic scintillators were placed in front of each detector to veto charged particles (such as 
decay electrons). The counters had a timing resolution of better than 5 ns full width half maximum (FWHM) for 
7-rays. The counters were about 30 cm from the target so a 7-ray time of fiight is 1 ns. 

Pulse shape discrimination was used to distinguish neutrons from 7-rays. Two different modules were used at 
different times, but the discrimination was always set conservatively Sources of ^"Co and AmBe were used to 
set up the system. If 7-rays are detected during the experiment, a prompt peak is clearly distinguishable, so it is 
straightforward to monitor the electronics on-line. 

The time of arrival of an event with respect to the stopping muon was measured by both a 5 /iS and a 10 /iS full 
scale time to digital converter (TDC). These TDCs have been tested and do not contribute significantly to the errors 
in the experiment. 

Data acquisition was carried out by a VAX station 3200 and a PDP-11 frontend processor (starburst). Data col- 
lected from the CAM AC modules were written to 8 mm tapes for later analysis. Further technical details of the 



2 



equipment are available in the thesis by Stocki [hd 



To test the overall system several targets were chosen for which there is no hyperfine transition, or for which the 
transition is too fast to be observed. The time of arrival of the neutrons could then be fitted by the formula: 

where D is the total disappearence rate for the clement under study, B is a flat background, to is the mid-point of 
the rise-time curve, and S is a folding of the time of flight spread and instrumental time resolution. 

A typical fit for Mg is illustrated in Figure |l} The values for S and to are presented in Table || and show the 
tendency for heavier elements to have a slower time of flight for the neutrons. As to refers to the time of arrival 
of the 7-rays, one can add the 7-ray time of flight (0.97 ns) to obtain the average neutron time of flight (At) and 
average neutron energy (E). From S one can remove the counter timing resolution to obtain the spread in energy of 
the neutrons (AE). These values were used as guidance in analyzing the elements with a hyperfine effect. Note that 
these values for E and AE are for the neutrons detected in this experiment, and do not constitute measurements of 
these parameters for muon capture in general. 

By taking the time derivative of the data such as that shown in Figure |], one can get an idea of the neutron time 
of flight and detector resolution effects. Figure |^ shows this time derivative spectra for Au, brass. Si, and S. One can 
clearly see in brass and Au that a single error function (or gaussian) is not enough. We found that the function: 

N(t)^ -erf^^ + - + R[ -erf^^ + -] 

seems to give a better fit for those elements. The effect is not seen in Si and S, due to limited statistics. 

For most elements there was a small correction (a few percent) for background from muons stopping in carbon in 
the scintillators. In addition, a few targets contained two elements viz LiF, (CF2)„, LiCl, and CCI4. We used teflon, 
i.e. (CF2)n, and CCI4 as test cases to ensure that the corrections could be adequately applied. For such molecular 
targets one needs to know the atomic capture ratio. The Fermi- Teller Z law is far from adequate and a Za variation 
gives a better approximation pT[ |. However, far more satisfactory is to use an actual measurement as the empirical 
variation is significant. 

For LiF the atomic capture ratio (Li/F) has been measured to be 0.28(3) by Zinov et al. jl^ and 0.10(8) by Wilhelm 
et al. |l3| for muons, and 0.10(1) Q and 0.22(2) for pions, which should be the same. As Zinov et al. disagree 
with later measurements for several other molecules, we take the pion values and use 0.16(12) with a liberal error. 
This value agrees with the empirical capture probability for muons determined by von Egidy and Hartmann which is 
0.18(5) for LiF For (CF2)„, Martoff et al. (15) found that (89.4 ± 0.7) % of pions capture on fluorine (i.e. an 
atomic capture ratio C/F of 0.24(2) in the normal definition, which takes into account the relative number of atoms. 
(The Z3 law gives 0.87; von Egidy and Hartmann |l^ ] do not give a value for carbon, but using the Martoff value, and 
a measurement of CO2 [|l^, we estimate the C capture probability is 0.33(9) in the von Egidy-Hartmann scheme). 
For CCI4 there is no measurement but using our value of the C capture probability we obtain C/Cl — 0.25(7) i.e. (5.8 
± 1.6) % of muons stop in the carbon of CCI4. The Zs rule gives C/Cl = 0.71, i.e. 15 % of muons stop in the carbon 
of CCI4. For LiCl the atomic capture ratio, as measured by Daniel et al. [|l^ is 0.19(8). 

Although these values are very uncertain, the situation is not as bad as it might appear. Since the lighter elements 
produce many fewer neutrons, the corrections turn out to be fairly minor. Hence, the number of stops in an element 
for which capture occurs (most of which produce neutrons), are 0.53% for Li, 7.8% for C, 33 % for F and 75 % for 
CI 1^. The effective background, therefore, from lithium in LiF is 0.25 %, which is negligible in comparison to other 
problems. We assume the neutron multiplicity to be the same for these elements. 

III. EXPERIMENTAL RESULTS 

The results for elements with a hyperfine effect are fitted to two different functions; one with a single error function 
viz 
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where H is the hyperfine rate, A is the hyperfine asymmetry, F and G are contributions from carbon in the counters, 
P and Q are contributions from other elements (if appropriate), and B is a flat background. The parameters F, G, P, 
and Q were calculated and fixed. These components make no substantial contribution to the data analysis. By using 
these two equations we can get a handle on the systematic errors. In Equation 4 we tie the Si and ti values to the S 
and to for two cases, using fits to the brass and the Au data. For LiF the data were of sufficient quality that all the 
other paramters (to, S, C, D, A, H) could be fitted freely. In this case the hyperfine effect is well separated from the 
rise-time effect caused by the neutron time of flight effects. The best fit is illustrated in Figure ^ and the values of S 
and to are given in Table [II. They are consistent with the non- hyperfine elements but slightly different. 



This is a salutary warning because we shall need to use both Equations (3) and (4) to get a handle on the systematic 
errors due to the time of flight and timing resolution effects. Our recommended values are the averages of these two 
methods. One would think that one could get the values of S and to from the time spectra of the neighbouring O"*" 
nuclei. This turns out to be a bad assumption. First let us understand the physics behind this hypothesis. Nuclei 
without a hyperflne effect are 0+ nuclei, often even-even nuclei such as ^^O, ^"^Mg, ^^Si, and ^°Ca. These are tightly 
bound. However odd Z nuclei like ^^F, ^'^Na, ^^Al, and '^^P tend to have N=Z-|-I and so, transforming a proton to 
a neutron proceeds to a nucleus even more neutron rich and more likely to fall apart. The typical neutron energy 
spectrum after muon capture is composed of two components, an evaporation spectrum peaked at about 1.5 MeV 
followed by a high energy tail starting around 5 or 6 MeV |l9j. Setting a threshold at 10 MeV Kozlowski et al. |2^] 
have shown that high energy neutrons constitute the following fraction: 26(5)% in ^^O, 19(3)% in ^*Si, 11(2)% in Ca 
and 11(2)% in Pb. Thus for light elements the high energy component is very important, and it is not surprising if 
it is sensitive to details of the nuclear structure of the product nucleus. The results of our experiment confirm that 
there are significant variations from nucleus to nucleus. 

For Na, Al, P, CI, and K we do not have sufficient data to allow a free parameter search, so we must average the 
single error function results with the two error function results, otherwise we do not know whether the results are 
from both time of flight and hyperflne effects. Notice for LiF and (CF2)„ this does not have to be done, because 
in F the capture rate asymmetry is large enough and the hyperflne rate is slow enough that the hyperfine rate can 
be distinguished from the time of flight effects. One way out of the dilemma for the other elements is to raise the 
energy threshold on the neutron detector. However, the hyperflne asymmetries are very small, and statistics becomes 
a problem. The other way is to do a full study of the neutron spectrum but this is a major project. We have therefore 
been forced to simply compare results from different fitting functions, using shape information from the spectra in 
Figure 2. The difference between the fits gives us an estimate of the systematic error. In Table [II we present some 
of these fits to illustrate the problem; Figure ^ shows two fits for Na which are both quite satisfactory. In Table IV 
our final recommended values are compared with previous results and the hyperfine asymmetries are compared in 
Table to various theories, and previous experiments. Note that in Table DVl the CI and K results did not use these 
error functions; since their asymmetry values were < 0, the hyperfine effect could be easily distinguished from the 
neutron time of ffight effects. The most important feature is that apart from ^^F and to some extent ^'^Na, the hyper- 
flne asymmetries are small, and that creates part of the problem. The 7-ray experiments observe large asymmetries 
(and 7-rays all travel at the same speed). Thus 7-ray experiments with lower statistics can effectively compete with 
neutron-detection experiments. A word of caution should be given concerning the comparison of the asymmetry with 
previous calculations, which were made for the total capture rates (including transitions to bound levels which give 
off no neutrons, and to levels which give off 2 neutrons). The experimental asymmetries are for speciflc experimental 
conditions with a threshold on the neutron energy and an energy dependent efficiency for the neutron detector. Thus 
one would not expect an exact equivalence, but the comparison is interesting. 

We see that, on the whole, the errors on the hyperflne rate are somewhat large. The two results which are most 
convincing are F and Na, but the latter is not sufRcient to resolve the disagreement between the earlier sodium metal, 
and NaF results. 
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The phosphorus result is worth a short discussion. An initial polarization measurement by Egorov et al. [ pT[ claimed 
that an asymmetry had been observed. The measurement was from to 2.4 /is implying no fast depolarization. Note 
also that the upper level is F=l and the lower level F=0, so that if a fast hyperfine transition does occur, there can be 
no residual polarization. Later asymmetry measurements by Hutchinson et al. and Babaev et al. observed no 
asymmetry and set limits at 10% of the value observed by Egorov et al., thus indicating that the hyperfine transition 
was fast, as suggested by Winston |^ who calculated 58 /is~^, i.e. r = 17 ns. This scenario was consistent with the 
data of Gorringe et al. |^ using 7-ray detection, but no positive observation was obtained. Again, our own result has 
systematic difficulties but is the first positive identification of the hyperfine transition in phosphorus. All the other 
hyperfine rates are consistent with earlier results. 

The comparision of Na and NaH in our data is interesting, but not of sufficient quality to resolve the questions 
about the sodium hyperfine rate, similarly for Al and LiAlH4, though the trend is for the hyperfine rates to be similar. 

In conclusion, we have measured the hyperfine rate in P for the first time, and obtained a more accurate value for 
F. Even though neutron detection appears to be the most sensitive technique for measuring the hyperfine transition 
rate, the problems with the neutron spectrum and the spread of the time of flight to the detectors have made it 
very difficult to observe fast rates. Further experiments should use a higher energy threshold (and take considerably 
greater statistics) . An alternative technique would be to use a 7-ray detector with a high resolution, but a faster time 
response than a HPGe detector. 
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FIG. 1. The time spectrum of neutrons from muon capture in Mg, fitted to equation 1. The data are presented with the 
muon disappearance rate divided out. 
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FIG. 2. The time derivative spectra of neutrons from muon capture in Au, brgiss, Si and S fitted to Equation 1 and Equation 
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FIG. 3. The time spectrum of LiF, with the decay divided out. The hyperfine transition in fluorine is clearly visible with its 
time constant of about 180 ns. 
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FIG. 4. The time spectrum of Na, with the decay divided out for the one error function case (top) and the two error function 
case (bottom). 
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TABLE I. 


Properties of the targets used. 






z 




Target 




Thickness 




Number of Events 










(g/cm^) 




in Time Spectrum 


3,9 




LiF 




2.70 




1.2 X lO"" 


6,9 




(CF2)„ 




2.25 




3.5 X 10^ 


11 




NaH 




1.14 




1.5 X 10® 


11 




Na 




1.65 




1.5 X 10'' 


12 




Mg 




1.61 




2.2 X lO'' 


13 




Al 




4.69 




1.2 X 10'' 


13 




LiAlH4 




1.51 




1.9 X 10^ 


14 




Si 




3.55 




3.0 X 10^ 


io 








2.79 




2.8 X iU 


17 




LiCl 




2.60 




7.6 X 10^ 


17 




ecu 




6.18 




4.9 X 10^ 


19 




K 




1.78 




8.4 X 10^ 


29,30 




Brass (Cu+Zn) 




2.79 




1.8 X 10^ 


TABLE II. 


Values of S and to from non-hyperfine targets. Also given are the average energy E, and the energy spread AE 


of the neutron events used for the timing measurements. 








Z 


Target 


S 


to 


At 


E (MeV) 


A E 






(ns) 


(ns) 


(ns) 


(McV) 


(McV) 


12 


Mg 


4.63 ± 0.24 


4.22 ± 0.23 


5.19 ± 0.23 


16.5 ± 1.6 


27 ±4 


14 


Si 


5.26 ± 0.17 


5.14 ± 0.17 


6.11 ± 0.17 


11.9 ± 0.7 


20 ±2 


16 


S 


4.84 ± 0.12 


4.35 ± 0.14 


5.32 ± 0.14 


15.7 ± 0.9 


25 ±2 


29,30 


Brass 


5.92 ± 0.11 


9.72 ± 0.11 


10.69 ± 0.11 


3.9 ± 0.1 


14.9 ± 0.7 




(Cu+Zn) 












79 


Au 


6.93 ± 0.08 


13.36 ± 0.07 


14.33 ± 0.07 


2.2 ± 0.1 


10.4 ± 0.3 
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TABLE III. Values of S and to for hyperfine targets, and their effect on the systematic uncertainties of If and A. The different 
results for H and A for each target should be taken as a measure of the systematic uncertainties. 



z 


Target 


S 
(ns) 




to 

(ns) 




Type of 

Error 
Functions 


(/^ 


H 

s-i 


) 


A 

(unitless) 


9 


LiF 


5.38 ± 0.18 
4.92 ± 0.15 
O.ol It u.zu 


(varied) 
(varied) 
(varied) 


5.30 ± 0.17 
4.81 ± 0.15 

0.0/ zt u.lo 


(varied) 
(varied) 
(varied) 






1 
2 
2 




5.6 ± 

5.7 ± 
5.6 ± 


0.2 
0.3 


0.288 ± 0.006 
0.290 ± 0.006 
0.287 ± 0.006 


9 


(CFa),. 


5.03 ± 0.31 
4.68 ± 0.26 

4.oy zt U.OD 


(varied) 
(varied) 
(varied) 


5.18 ± 0.29 
4.59 ± 0.27 

0. / 1 zt U.Ol 


(varied) 
(varied) 
(varied) 






1 
2 
2 




5.2 ± 
5.2 ± 
5.2 ± 


0.4 
0.4 
0.4 


0.33 
0.33 
0.33 


± 0.01 
± 0.01 
± 0.01 


11 


NaH 


5.0 ± 0.2 ( 

4.31 ± 0.22 

A QA -1- n 1 

4.y4 zt U.ID 


varied) 

(varied) 

(varied) 


5.0 ± 0.1 ( 

4.51 ± 0.13 


varied) 

(varied) 

(varied) 






1 
2 
2 




26 
32 
26 


± 
± 
± 


4 
4 
5 


0.15 
0.19 
0.15 


± 0.02 
± 0.03 
± 0.02 


11 


Na 


4.7 ± 0.1 ( 

A 99 -i- n 1 Q 

4.ZO zt U.lo 

5.21 ± 0.23 


varied) 

(varied) 

(varied) 


4.5 ± 0.2 ( 

9 77 -1- n '^0 
0. t i zt U.oZ 

4.75 ± 0.14 


varied) 

(varied) 

(varied) 






1 
2 
2 




oy 

49 
34 


± 7 
± 10 
± 6 


0.18 
0.24 
0.15 


± 0.04 
± 0.06 
± 0.03 


13 


Al 


5.13 ± 0.13 
4.63 ± 0.13 
5.21 ± 0.23 


(varied) 
(varied) 
(varied) 


3.52 ± 0.38 
2.73 ± 0.30 
4.75 ± 0.14 


(varied) 
(varied) 
(varied) 






1 
2 
2 




68 
75 
71 


± 
± 
± 


9 
7 
9 


0.30 
0.37 
0.31 


± 0.06 
± 0.04 
± 0.06 


13 


LiAlHi 


5.04 ± 0.31 
4.96 ± 0.23 


(varied) 
(varied) 


2.37 ± 0.45 
4.66 ± 0.85 


(varied) 
(varied) 






1 
2 




90 
53 


± 6 
± 25 


0.50 
0.22 


± 0.04 
± 0.13 


1 


p 


5.90 ± 0.07 
6.73 ± 0.11 


(varied) 
(varied) 
(varied) 


5.7 ± 0.23 
4.57 ± 0.19 
5.45 ± 0.20 


(varied) 
(varied) 
(varied) 






1 
2 
2 




44 
62 
76 


± 
± 
± 


8 
6 
5 


0.12 
0.23 
0.27 


± 0.03 
± 0.03 
± 0.02 








TABLE IV. Muonic hyperfine transition rates 


in ^s ^. 












Compound 




Theory | 






Previous Work 








This experiment 


LiF 

NaF 

(CF2)„ 






5.8 






5.8 ± 0.8 
6.3 ± 1.8 [ 

4.9 ± 1.2 [ 


2 
I 


5.6 ± 0.2 
5.2 ± 0.4 


Na 
NaH 
NaF 






14 






15. 5± 1.1 
8.4± 1.9 


1 


1 










38 ±9 
28 ±5 


Al 
LiAlH4 






41 






41 ± 9 1 














72 ±9 
89 ± 10 


Red P 






58 






Ah > A_ 
A- > Xh 














60 ± 15 


LiCl 
CCI4 






8 






6.5 ± 0.9 














14 ± 27 

15 ± 29 


K 






22 




















25 ± 15 



''A- f» 1.1 /is-i 
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TABLE V. The capture rate asymmetry. The theory values are for aU captures. This experiment is for neutron detection 
only (weighted by multiplicity). 



Compound 




r^rimaKOit [L | 


T T"l--v ^ 1 1 

U Derail 


25|| 


Previous work 1 2 ] 


This Experiment 


LiF 


0.24 


0.36 


0.36 




0.36 ± 0.04 (n) 

U.zo ± U.U4 ('y ) 


0.29 


± 0.01 


(CF2)„ 


0.24 


0.36 


0.36 






0.33 


± 0.01 


i\a 


U.Uo 










0.17 


± 0.05 




U.Uo 










0.16 


± 0.03 


Al 


0.09 


0.14 


0.22 






0.34 


± O.Oo 


LiAlH4 


0.09 


0.14 


0.22 






0.35 


± 0.15 


Red P 


0.16 


0.22 


0.25 






0.2 


± 0.1 


LiCl 


-0.06 ('''■CI) 
-0.06 (^■'Cl) 


-0.10 (•''CI) 








-0.03 


± 0.06 


ecu 












-0.03 


± 0.06 


K 


-0.05 


-0.07 








-0.080 


± 0.055 
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